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Abstract. We propose that radio-quiet quasars and Seyfert galaxies have central black holes powering outflows and jets which 
propagate only for a short distance, because the velocity of the ejected material is smaller than the escape velocity. We call them 
"aborted" jets. If the central engine works intermittently, blobs of material may be produced, which can reach a maximum radial 
distance and then fall back, colliding with the blobs produced later and still moving outwards. These collisions dissipate the 
bulk kinetic energy of the blobs by heating the plasma, and can be responsible (entirely or at least in part) for the generation 
of the high energy emission in radio-quiet objects. This is alternative to the more conventional scenario in which the X-ray 
spectrum of radio-quiet sources originates in a hot (and possibly patchy) corona above the accretion disk. In the latter case the 
ultimate source of energy of the emission of both the disk and the corona is accretion. Here we instead propose that the high 
energy emission is powered also by the extraction of the rotational energy of the black hole (and possibly of the disk). By means 
of Montecarlo simulations we calculate the time dependent spectra and light curves, and discuss their relevance to the X-ray 
spectra in radio-quiet AGNs and galactic black hole sources. In particular, we show that time variability and spectra are similar 
to those observed in Narrow Line Seyfert 1 galaxies. 
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1. Introduction 

One of the most popular scenario to explain the dichotomy 
between radio-loud and radio-quiet Active Galactic Nuclei 
(AGNs) assumes that only rapidly spinning black holes can 
give rise to the relativistic jets responsible for the radio emis- 
sion and higher frequency non-thermal radiation observed in 
radio-loud objects (e.g. Blandford 1990). Therefore it was with 
some surprise that the first evidence for rapid rotation of a 
black hole came from a radio-quiet object, namely MGC-6- 
30-15. In the X-ray spectrum of this Seyfert 1 a broad iron line 
was observed by ASCA and BeppoSAX (Tanaka et al. 1995; 
Guainazzi et al. 1999) to be consistent with emission outside 
6i? g (R g = GM/c 2 ), i.e. the innermost stable orbit of an ac- 
cretion disc around a static black hole. Recent XMM-Newton 
observations (Wilms et al. 2001; Fabian et al. 2002), however, 
indicate that the emission may extend well within 6i? g , thus re- 
quiring a rotating black hole (the innermost stable orbit in the 
case of a maximally rotating black hole being R g ), and con- 
firming the result obtained by Iwasawa et al. (1996) during a 
low flux state of the source observed with ASCA. Moreover, 
the very steep radial dependence of the iron line emissivity was 
interpreted by Wilms et al. (2001) as evidence for the extraction 
of the spin energy of a Kerr black hole, even if a pure geomet- 
rical explanation (but still requiring a rotating black hole) is 
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possible (Martocchia et al. 2002), provided that the illuminat- 
ing source is very close to the black hole and resides on the 
symmetry axis. Similar evidence comes from XMM-Newton 
observations of the galactic black hole candidate XTE J 1650- 
500 (Miller et al. 2002a). 

Furthermore, Elvis, Risaliti & Zamorani (2002) have re- 
cently suggested that the X-ray background requires a high ef- 
ficiency of mass to energy conversion in the accretion process, 
possible if the black hole is rotating, but problematic in the case 
of a Schwarzschild black hole. 

These results are at odds with the idea of a slowly spinning 
black hole in radio-quiet objects. On the other hand, and in a 
complementary way, radio-quiet objects are not radio-silent: 
even if the "dichotomy" between radio-loud and radio-quiet 
objects is currently under scrutiny (see e.g. White et al. 2000 
and Ivezic et al. 2002), all AGNs can produce radio emission at 
some level, which in turn is consistent with the idea that some 
sort of jet or outflow is always present, responsible to acceler- 
ate electrons to relativistic energy to radiate by the synchrotron 
process in the radio band. This idea has received recently fully 
support by VLBI imaging of Seyfert galaxies (Ulvestad 2003, 
and references therein) which revealed the presence of a mini- 
jet (at the sub-pc scale) in many Seyfert (radio-quiet) galaxies. 
In several cases it was also possible to detect the proper mo- 
tions of knots in the jet, which appear to move with subluminal 
apparent velocities of the order of a tenth of the speed of light. 
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Therefore it is conceivable to assume that all black hole 
plus accretion disk systems in AGNs can produce some kind 
of outflow or jet, but that only in a minority of cases (i.e. the 
"pure" radio-loud objects) the jet is successfully launched and 
accelerated to relativistic speeds. In the majority of cases, the 
jet is "aborted", yet it is responsible for a relatively weak ra- 
dio emission. The idea that all AGNs produce a jet is not new: 
among others, Falcke & Biermann (1995) suggested the jet- 
disk symbiosis for all AGNs, while Henri & Petrucci (1997) 
and Malzac et al. (1998) have argued that the initial part of a 
jet in a radio-quiet object can produce relativistic particles illu- 
minating the disk. In these scenarios, however, the jet has either 
bulk relativistic motion or it contains very energetic particles, 
which are relativistic in the jet-comoving frame. In our sce- 
nario, as it will be explained below, the jet has sub-relativistic 
bulk velocities by assumption, and also most of the emitting 
electrons are thermal, with subrelativistic temperatures. 

It is possible that a source which is usually radio-quiet may 
occasionally be successful in launching relativistic jets. This 
could explain the properties of galactic superluminal sources, 
in which major outflows sometimes occur. If this is true, these 
sources should be considered a crucial link between radio-loud 
and radio-quiet objects. In this respect, it is worth noting that 
in Galactic superluminal sources and Galactic black hole can- 
didates there is often (even in radio-quiet states) the presence 
of a high energy X-ray power law, which may be associated 
with the emission from a jet, or at least from an outflow (as in 
the case of XTE 1 1 18+480: Miller et al. 2002b). 

The aim of this paper is twofold. First we will explore if 
the simplest "abortion mechanism" which comes into mind, 

1. e. a "jet" which does not succeed to reach the escape velocity, 
can work, at least qualitatively. Then we explore the possibility 
that the power initially in the jet and/or outflow can be used 
to heat the particles responsible to emit the X-ray flux from 
radio-quiet AGNs. In other words, we substitute the popular 
hot corona, possibly patchy, which sandwiches the accretion 
disk (i.e. Haardt & Maraschi 1991), with a single hot region on 
the rotation axis of the black hole, thought to be the site of the 
jet abortion. We then perform numerical simulations assuming 
to launch many blobs with slightly different velocities and time 
separations, calculate their trajectories and follow their evo- 
lution, accounting for the collisions occurring between them. 
This allow us to calculate the produced luminosity in each col- 
lision, and the total luminosity received by the observer in the 
likely case that more than one shell-shell collision is occurring 
at any given observing time. We will also study if the typi- 
cal scattering optical depths and temperatures of the scattering 
particles in the aborted jet scenario are in agreement with what 
observed (Petrucci et al. 2001; Perola et al. 2002). Finally we 
will discuss our findings and derive some observational conse- 
quences enabling to test this scenario. 

2. Escape velocity from a Kerr black hole 

The mechanism to form, accelerate and collimate jets in radio- 
loud objects is not well understood, even if several propos- 
als exist in the literature (see e.g. the review by Lynden-Bell 
2001). Here we do not attempt to propose a new mechanism. 



We simply postulate that in radio loud-quiet objects a similar 
mechanism is at work, but it is, on average, not able to impart 
to the outflowing material a bulk velocity larger than the escape 
velocity. We consider a blob of material in ballistic motion, and 
neglect the influence of magnetic fields and/or radiative fields 
(for calculations including the effect of accretion disk radiation 
see Vokrouhlicky & Karas 1991). 

The equation of motion of a test particle along the rotation 
axis of a Kerr hole is given by (see e.g. Vokrouhlicky & Karas 
1991): 



dx 



f3j 2 (x 2 
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where x = R/R g , a is the (dimensionless) specific angular 
momentum (~ 1 for a maximally rotating Kerr hole), and R g = 
GM/c 2 is the gravitational radius. Its solution is given by: 
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It is then possible to find the escape velocity for a test particle 
having an initial velocity c/3q at the distance xo = Rq /R g on 
the rotational axis of a Kerr hole: 
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yielding /3 CSC = 1 for R = R g and a = 1 (maximal Kerr 
hole), and for R = 2R g and a = (Schwarzschild hole). Fig.^ 
shows the escape velocity /3 CSC as a function of the initial radius 
xq for different values of the angular momentum parameter a. 
Material with (3 < (3 CS c will then outflow in the radial direction 
reaching a maximum distance from the black hole, and then it 
will fall back. The maximum distance above the hole where the 
motion inverts, x max , can be found setting 7(x max ) = 1 in Eq. 
(fj), to find 



k+ y/k 2 (l-a 2 ) + 2ak 
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xl + a 2 



(4) 



-ll{xl +a 2 - 2xo) 
Fig.|2]shows for illustration some examples of "trajectories". 

3. Jet power and internal shocks 

In this section we discuss the basic features of our model. In or- 
der to have an analytical treatment, we will consider here only 
a simple case, while the more realistic, but much more complex 
case, will be treated numerically in section 5. 
The total power carried by the jet is 



Lj = r M out c 2 



(5) 



where To is the bulk Lorentz factor corresponding to the initial 
velocity of the jet, and M Q ut is me mass ejection rate, which is 
in principle a free parameter of the model, but it can be con- 
strained by the observed ratio between the X-ray and UV lu- 
minosity. A fraction ijj of Lj will be transformed into radiation. 
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Fig. 1. Escape velocity as a function of distance form the black 
hole and for a =0, 0.7, 0.9, and 1, from top to bottom. 
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Fig. 2. Trajectories of test particles in motion along the ro- 
tation axis of an extreme Kerr hole (a = 1, solid lines), 
and Schwarzschild (dashed lines), for the initial radii xq = 
R/R g =5, 5.5, and 6, as labeled. The initial velocity is /3o = 0.5. 
Time is measured in units of the initial light crossing time 
Ro/c. 



We propose that the central engine producing the jet is not 
working continuously, but intermittently, leading to collisions 
between blobs of material launched at different times, and mov- 
ing with slightly different velocities. In the case of success- 
ful jets this will lead to the formation of "internal shocks", 
first proposed by Rees (1978) for the jet of the radio-galaxy 
M87, later proposed to explain the prompt emission of gamma 
ray bursts (Rees & Meszaros 1992; Rees & Meszaios 1994; 
Sari & Piran 1997), and recently proposed again to explain 
the entire spectral energy distribution of all radio-loud AGNs 
(Ghisellini 1999; Spada et al. 2001). In this scenario the shells 



collide at a typical distance i?diss from the black hole given 
by i?diss ~ r 2 i?o> where T is the average bulk Lorentz factor 
of the shells and i?o is the initial shell-shell separation. In the 
case of aborted jets the intermittency of the process will lead to 
collisions between a blob already falling backwards to the hole 
and a blob still moving upwards. There is no need, in this case, 
to assume slightly different velocities, since the collisions will 
occur anyway. The details of the collision depend upon the ini- 
tial velocities and the typical "launching" site, which is related 
to the time needed to accelerate a blob. 

There may be a range of locations i? co ii where collisions 
preferentially take place, depending on the initial time separa- 
tion of two consecutive shells, their initial velocity /?o, and on 
the initial launching site Rq. These parameters, together with 
the angular momentum of the Kerr hole, determine R co \\, which 
in turn determines the amount of bulk kinetic energy which can 
be dissipated. Larger R co \\ in fact correspond to smaller kinetic 
energies (and larger potential ones). Fig.|5]shows for illustra- 
tion simple examples of the velocity and kinetic energy of two 
colliding shells, as a function of i? co ii- In these examples we 
assume the same mass, initial velocity and launching site for 
all blobs. 

In a more realistic case, R co \\ will be dependent partly on 
the exact values of the initial parameters, and partly on the past 
history of the process, in the sense that a falling blob may inter- 
act with more than one later blobs (including some which may 
have (3 > f3 csc ). Interestingly, the power which can be extracted 
is a function of the dissipation distance, since the more power- 
ful collisions will be between just launched blobs and falling 
blobs with a back velocity of the order of the initial one, and 
this will occur closest to the hole (see Fig. |3j. We also stress 
that in these examples we have only treated the simplest possi- 
ble case, neglecting in particular the interaction of the blob with 
the radiation produced by the accretion disk and the motion off 
the angular momentum axis. 

With our simple approximations, the dissipated power is 
the entire kinetic power of the shells when they collide, since 
they have equal and oppositely directed velocities. Thus the 
total momentum is zero if the shells have, as we assume, equal 
masses. We can then link the radiative efficiency of the process 
to the kinetics of the collision through: 



where e c is the fraction of the collisional energy radiated by 
the electrons (and possibly electron-positron pairs), and r co n 
is the bulk Lorentz factor when the blobs collide. We do not 
know what is the dominant acceleration mechanism transform- 
ing ordered into random energy, and this precludes to know 
with confidence what fraction of the total energy goes to the 
emitting electrons (and possibly positrons), and what fraction 
is instead given to protons and to magnetic field. Equipartition 
between protons, electrons and magnetic field would result in 
a value of e c ~ 1/3. The value of r/j can be calculated for any 
assumed value of e c from Fig.|3]and Eq. (|6}. 

Observationally, the X-ray luminosity is of the order of 10- 
50 per cent of the luminosity in the optical-UV (e.g. Walter & 
Fink 1993), thought to be produced by the accretion disk, for 
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which Ld = ?7dMi„c 2 . Here M- m is the mass accretion rate and 
r^d is the accretion efficiency. 

If all the X-ray luminosity Lx comes from the aborted jet, 
i.e. Lx = VjLj, while ~ £tjv> we have: 



M 



1 Vd L x 



(7) 



Since T is always of order unity, we have that M out ~ 
(Lx / L\jy)Mi n ~ 0.1Mi n if the jet and the accretion efficien- 
cies are of the same order. 



4. The X-ray emission 

In this section we again discuss the very simplified case of the 
collision of two oppositely directed blobs having the same mass 
and same velocity (in modulus), in order to understand the ba- 
sics of the interaction with an analytical treatment. 

4. 1. Formation of the X-ray continuum 

Having established that the simple aborted jet scenario dis- 
cussed above can in principle account for the energetics of the 
X-ray emission, we must now check whether the physical pa- 
rameters of the aborted jets, i.e. its optical depth and tempera- 
ture, are consistent with the X-ray spectral constraints. 

The accelerated leptons are embedded in a very large radi- 
ation energy density produced by the close-by accretion disk 
and by the local magnetic field possibly enhanced by the col- 
lision. Under these conditions, the most efficient radiation pro- 
cesses are Inverse Compton scattering and cyclo-synchrotron 
emission. We will see below that energy balance ensures that 
the equilibrium lepton energy is mildly relativistic at most. This 
implies that the cyclo-synchrotron process occurs in the self- 
absorbed regime and then does not contribute to the cooling. 
This process could however be important for establishing an 
electron Maxwellian distribution in a timescale shorter than 
the dynamical time (Ghisellini, Haardt & Svensson, 1998). For 
simplicity, we then assume that the energized leptons have a 
thermal energy distribution. It is worth noting that a distribu- 
tion which is not a perfect Maxwellian, but has a well defined 
mean energy, such as a distribution with a well defined peak en- 
ergy, or a power law energy distribution (oc E~ n ) with n > 3, 
give rise to a Comptonization spectrum which cannot be dis- 
tinguished by the one formed by a perfect Maxwellian (see e.g. 
Ghisellini, Haardt & Fabian 1993). 

To find the typical average energy 7m c 2 of the emitting 
leptons, we balance radiation losses with the energy gains due 
to the shell-shell collision process. The inverse Compton cool- 
ing rate is: 



E r . 



-a T cU vl 2 (3 2 
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where U T is the radiation energy density of the disk emission. 
The heating rate due to the shell-shell collision is 
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Fig. 3. Consecutive shells collide at different distances from 
the black hole and with different velocities. In the collision they 
dissipate their kinetic energy, proportional to r co n — 1, An ini- 
tial distance xq = 5 and initial velocity (3q = 0.6 is assumed. 
Solid lines refer to a maximally rotating Kerr hole, dashed line 
to a Schwarzschild black hole. Top panel shows the velocity 
of the shells when they collide; the mid panel shows the total 
dissipated energy, which is greater closer ot the hole; the bot- 
tom panel shows the Comptonization y parameter, calculated 
following Eq. H2\ . and assuming, as labeled, Lj = £d and 
e c = 0.3. We show the case of shells of fixed dimensions (set 
to 5i? g in this example) or expanding as r 3 = ijjR co \\, where 
we have used ip = 30°. 

where N c is the total number of leptons in the emitting region. 
Balancing heating and cooling we have: 



7 2 /3 2 - 3 ^ 



4 OrcNJJr 



(10) 



Assume that the "jet" emitting region has a transverse size rj 
and a width Arj . Its scattering optical depth is 



frfAr: 



(11) 



where / is a geometry dependent factor (e.g., / = 47r/3 for a 
sphere). 

If tt<1, the Comptonization parameter is y = 
(4/3)tt7 2 /3 2 . In this case we have: 



(9) y 



frfcU r 



(12) 
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The energy density U T depends on the specific geometry of the 
system and on the accretion luminosity. Note that Ld should 
always be greater than rjjL j/2, since about half of the jet pro- 
duced luminosity impinges on the accretion disk. However, ob- 
servations typically give La ~ Ljjy ~ lOLx, suggesting that 
the accretion-produced luminosity is dominating with respect 
to the jet luminosity intercepted and reprocessed by the disk. 
For simplicity, let us assume a Newtonian disk producing a to- 
tal power id, which has a minimum radius r*. Per unit surface 
area the dissipated flux is (e.g., Frank, King & Raine 1985): 



D(r) 



3Ld^* 
4 7rr 3 



I T * X 1 2 
r 



(13) 



where r is the radial coordinate of the accretion disk. The inten- 
sity J(r) corresponding to this dissipation is I(r) = D(r)/%. 
At the location R co \i the total radiation energy density produced 
by the disk is 



C^r(-Rcoll) 



3L<ir*R, 



coll 



2ttc 



1 - (r./r) 



1/2 



(^oll + r 2 ) 



3/2 



dr 
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where /i = cos 8 and 8 is the angle with respect to the normal 
of the accretion disk. 

Inserting f7 r (i? co n) in Eq. d!2t we can calculate the 
Comptonization y parameter. In the bottom panel of Fig. [3] we 
show, as illustration, how y changes as a function of i? co n if 
the shell remains of the same dimension (i.e. rp constant) or 
if instead it expands as a cone of semiaperture angle ip (in this 
case we have assumed rj = tpR co \\). In both cases we show the 
behavior for a maximally rotating Kerr hole (a ~ 1) and for a 
Schwarzchild hole (a = 0). 

A constant r- 3 implies an almost constant y-parameter. This 
is because rjj and U r decrease with i? co ii in approximately the 
same way. If the shells expand, instead, y decreases because 
the optical depth of the scattering electrons decreases as tt oc 



rf 2 oc R: 2 



■-coir 



Eq. d!2i also gives the j/-parameter in the classical case 
of an homogeneous corona sandwiching the accretion disk. 
Assuming that the jet dissipates the entire available gravita- 
tional energy, and that half of it is reprocessed by the accre- 
tion disk, we can substitute i]jLj with Lx, id with Lx/2 
Furthermore, the disk radiation energy density in this case can 
be approximated by U T = 2Lx/ (irr 2 c). Setting / = ir/2, we 
finally obtain y = 1. 

4.2. Are electron-positron pairs important? 

The assumption that the jet carries a power ij in the form of 
kinetic energy links the amount of transported power with the 
scattering optical depth. On the other hand, tt may not be lin- 
early proportional to Lj, because of the possible presence of 
electron positron pairs, which contribute fully to the scattering 
optical depth but not so much to the jet power (if protons are 
also present and are dominating the jet inertia). We will there- 
fore estimate first a lower limit to the scattering optical depth 



assuming no pairs, and then we will make some considerations 
to evaluate the contribution of pairs. 

To calculate the initial optical depth of each blob, without 
the contribution of pairs, we first consider the kinetic energy 
carried by each blob, and use it to find its mass Mf. 



L-,At 

Mi = — 



(15) 



where At is the time interval between two consecutive blob 
ejection. All the values of the quantities in Eq.^]have to be 
considered as average values. The total number of electrons 
contained in each blob is therefore iV c = Mi/m p if the jet is 
made by an electron-proton plasma. Inserting this value in Eq. 
II llwe get: 



cr T LjAi 



4tt U 



m p c 2 fr 2 T ii0 c 2 



f^i,0 -^Edd 



At 

r s/ c 



(16) 



For the second equality, we have used the jet luminosity in units 
of the Eddington one, the blob size in units of r g and time in 
units of r-g/c. In this way it becomes clear that the optical depth 
(without the contribution of pairs) is scale invariant with re- 
spect to the mass of the black hole. Furthermore, it is also clear 
that, for rj of the order of the gravitational radius, for jet ki- 
netic powers of order of 0.01-0.1 of the Eddington luminosity 
and for time intervals of a few times r g /c the initial optical 
depth is of order of 0.01-0. 1 . 

If the jet, instead, is made by a pure electron-positron 
plasma, then the initial value of the optical depth is a factor 
m p /m e greater, and consequently tt becomes much larger 
than unity. In this case the pairs annihilate efficiently, since the 
annihilation timescale is of order of ^/(ctt), and most of the 
kinetic power of the jet is lost through annihilation. We con- 
clude that if the jet is energetically important, then its inertia 
must be given mainly by protons. 

We stress that there is a remarkable feature in the model, 
which is the link between the kinetic power and the amount of 
electrons and protons carried by the jet. In pair-corona mod- 
els, the optical depth of the corona has a lower limit given by 
the created pairs (Haardt & Maraschi 1991). Here, instead, the 
limit is given by the jet power. 

We can also estimate the possible importance of pairs cre- 
ated during the emission phase. In other words, when primary 
electrons are heated by the shell-shell collision, they can emit 
photons above the pair production threshold: if the source is 
sufficiently compact, these photons create pairs which increase 
the optical depth and contribute to the emission. 

The effects of pairs have been studied in detail assuming a 
steady source in pair equilibrium (creation equal annihilation, 
Svensson 1984). Our source is probably never in steady state, 
but we can use the results of these studies as a guide to estimate 
the relevant pair optical depths and temperatures. The optical 
depth due to pairs, produced in steady state and in pair equi- 
librium without pair escape, can be approximated by (Haardt 
1994): 



T± 



0.016 (4) 



0.56. 



20<4<200 



(17) 
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where £ K = (JtL^/ '(rjm e c 3 ) is the compactness of the X-ray 
radiation produced by a shell-shell collision, and can be con- 
sidered equal to r/jij. This equation assumes that the radiation 
process is Comptonization, and is valid assuming spectral in- 
dices a x of the Comptonized spectrum around unity. The range 
of X-ray compactnesses relevant for Seyfert galaxies is within 
the same range of the compactness relevant for the accretion 
disk emission, i.e. £ x ~1— 10. 

By comparing Eq.^]with Eq.[H)]we see that pairs should 
be unimportant for compact jets. Furthermore, as will be clear 
in the following section, the optical depth of each shell is in- 
creasing each time it collides with another one (if lateral ex- 
pansion can be negleted), driving the typical optical depth to 
values a factor ~10 larger than the initial value. We then con- 
clude that electron positron pairs do not play a fundamental 
role in this scenario. 

5. Numerical simulations 

We have discussed so far the illustrative and very simple case 
of a pair of blobs of equal mass, equal launching site and initial 
velocity, colliding at different distances from the black hole. In 
a more realistic case, we should assume that the time interval 
At between the launch of consecutive blobs is variable, as well 
as the blob mass, initial velocity, and their launching site. To 
simulate this, we use a Monte Carlo code which extracts the 
initial quantities within an assigned distribution. We have then 
assumed a Gaussian distribution for the initial velocity (i.e. we 
assign a mean value and a width a), while for the time interval 
At we assume a Poisson distribution. We have then considered 
the case of a Gaussian distribution of the blob masses, leading 
to initially different kinetic energies of the blobs, and the case 
of equal kinetic energies for all blobs. In this latter case the 
mass is therefore determined by the blob velocity. 

We follow the trajectory of all blobs, and when they collide 
we calculate the corresponding dissipation through the equa- 
tion of conservation of energy and momentum: 

TiTOi + r 2 m 2 = Tf (mi + m 2 + e'/c 2 ) (18) 

Tiftmi + T 2 (3 2 m 2 = T f (3 f (mi + m 2 + e'/c 2 ) (19) 

In these equations the subscripts 1 and 2 stand for the two 
blobs, T = (1 — /3 2 )~ 1//2 is the bulk Lorentz factor just before 
the collision and V f if the final Lorentz factor of the two blobs 
just after the collision, e' is the dissipated energy, calculated in 
the frame comoving with the merged blobs. Hereinafter primed 
quantities are calculated in this frame. A completely anelastic 
collision is assumed. The two unknowns are j3f (or T f) and e' 
which are given by: 

TxP\mi + T 2 /3 2 m 2 

Pi ; 



Timi + T 2 m 2 



(20) 



where t' diss is the time needed to radiate the energy e e e'. In the 
dense seed photon environments we are considering, the elec- 
tron radiative cooling time is fast, always shorter than blob- 



blob crossing time t' c 



1/2 



^2 = [m 1 +m 2 + 2miTO 2 rir 2 (l - /3i/3 2 )J -mi-m 2 (21) 

A fraction e c of e' is given to the radiating electrons, making 
the X-ray luminosity, as observed in the comoving frame: 



cross- We then assume that t' diss 



tL 



We compute t' cl 

''diss — cross 



s in the frame of one of the two blobs, 

An 1-0102 



C \0X-02 



(23) 



We identify L' x with the "radiative jet luminosity" r?jLj 
of Eq. (I12> and derive the corresponding value of the 
Comptonization parameter y. 

At any time, the observer can see more than one collision, 
and therefore we calculate, at any time, the number of collid- 
ing shells. We then sum up their luminosities and assume, for 
each collision, a triangular luminosity profile with simmetric 
rise and decay timescales (each of duration idissX whose time 
integral is equal to the energy dissipated by the electrons. The 
observed luminosity produced by each collision is modified by 
special relativistic effects, which we account for introducing 
the Doppler beaming parameter: 

5 = iyi-ftcosfl) (24) 

where 9 is the viewing angle and fif must be considered with 
its sign (positive for shells approaching the observer). We then 
multiply the comoving luminosity by <5 4 and divide the intrinsic 
timescale by S. Note that S can be greater or smaller than unity, 
depending on 0f and 9. We take into account time dilation due 
to general relativity, by dividing intrinsic times by the factor 
(1 — 2-Rg/i?) 1 / 2 , where R is the distance from the black hole, 
and multiplying intrinsic frequency by the same factor. We call 
the resulting observed luminosity L° hs . 

The luminosity is assumed to be produced when two shells 
merge. For simplicity, we assume that the entire produced ra- 
diation is emitted once the velocity of the merged shell has al- 
ready reached the final value, and use this velocity to calculate 
the appropriate special relativistic effects. We then continue to 
follow the merged shells until a new collision occurs, or until 
it reaches a distance from the black hole equal or less than the 
initial one. 

Note that the mass of a generic blob increases for each col- 
lision, and in the absence of side expansion this implies a cor- 
responding increase of the scattering optical depth. 

We have also taken into account the additional electron 
heating due to Coulomb collisions between protons and elec- 
trons, occurring on the timescale t' ep given by (e.g. Stepney 
1983) 



, Ar^ m p /m c 

CTt 



In A 



m e cr 



3/2 



(25) 



where T c and T p are the electron and proton temperatures, and 
In A ~ 20 is the Coulomb logarithm. Coulomb collisions be- 
tween proton and electrons allows electrons to be heated for a 
time t cp . The additional energy gained by the electrons is radi- 
ated on a cooling timescale. This corresponds to an additional 
luminosity 
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100 200 300 400 500 
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Fig. 4. Light curves of the integrated X-ray luminosity (in units 
of the Eddington one) corresponding to numerical simulations, 
as discussed in the text. Thin lines corresponds to the power 
emitted by shell-shell collisions. To this, we have added a con- 
stant luminosity (0. 1 % Eddington) and the sum is shown by the 
thick solid lines. Note that these light curves take into account 
the total number of shells which an observe sees "on" at a given 
time. The three cases differ only for the average time interval 
between the blobs, which is 4, 8 and 16 Rq/c from top to bot- 
tom. Note that these results are scale-invariant: the light curve 
does not change changing the black hole mass. 

where e p is the fraction of dissipated energy heating the pro- 
tons. We then transform L' x p as before, taking into account 
beaming and gravitational redshift. In our cases t' CTOSS is almost 
always shorter than t' cp , and consequently L' x p is almost al- 
ways smaller than L' x , but it lasts longer, somewhat smoothing 
the lightcurve. 

For the simulations shown in Fig. and in Fig. [5] we have 
assumed: 

- (A)) = 0.5; <7/3 — 0.2; Gaussian distribution; 

- (At) =4, 8, 16 Rg/c; Poisson distribution; 

- launching site i?o = 5i? g equal for all blobs; 

- blob size rj = AR g equal for all blobs; 

- blob width Arj = 7-j equal for all blobs; 

- the blobs do not expand nor contract; 

- the jet initial kinetic luminosity Lj equal to the radiative 
accretion disk luminosity, both being equal to 0.1 ^Edd; 

- maximally spinning black hole; 

- e c = 1/3; e p = 1/3 for all blobs; 

- cos6> = 0.5. 

Clearly, some of these assumptions are not fully justified 
physically, but have been done just for ease of computation. In 
particular, one expects that the blob size and width change dur- 
ing the blob trajectory; one furthermore expects that the launch- 
ing size is not the same for all blobs. 




10" 5 0.0001 0.001 0.01 
log(L x /L Edd ) 

Fig. 5. Thomson scattering optical depth, plasma temperature 
(in keV) and Comptonization y parameter as a function of in- 
tegrated X-ray luminosity produced in single collisions for the 
case corresponding to the top panel of Fig.|4] The dashed verti- 
cal line corresponds to 0. 1 % of the Eddington luminosity which 
we have assumed to be constant. 



0.6 



0.4 



0.2 



~i — i — i — i — i i i 




Fig. 6. When a blob is launched, it has an optical depth be- 
longing to the distribution shown by the dotted histogram. The 
other histogram shows the distribution of tt soon after each 
collision. 



Fig- E] shows the light curves of the bolometric luminosity 
for three cases, corresponding to (At) = 4, 8, 16 R g /c (from 
top to bottom). Being an intrinsically intermittent process, the 
resulting light curves (thin lines) show a very large range of 
variability (the luminosity can virtually vanish for some time 
interval). This corresponds to the case of a completely "on and 
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Fig. 7. Light curves of the integrated X-ray luminosity (in 
units of the Eddington one) corresponding to numerical sim- 
ulations. The three cases differ for the average initial aver- 
age velocity, which is 0.3, 0.5 and 0.7 c, from top to bottom. 
The average time interval between consecutive shell ejection is 
(At) = 4i?Q- All other parameters are the same as in Fig. [4] 

off mechanism, i.e. there is a vanishing particle density be- 
tween the launched blobs. We consider this unphysical, even 
if useful for estimating the radiated power in the way we have 
discussed above. A more sophisticated treatment should take 
into account a "smoother" behavior of the central engine, with 
density and velocity profiles described by smooth functions. 
While this is referred to future work, the general effect on the 
light curve will be to allow some of the energy carried by the 
blobs (that in a more physical scenario can be thought as over- 
dense regions) to be dissipated all along the jet, and not only 
during collisions with other overdense region. This emission 
may correspond to a minor fraction of the luminosity produced 
during each collision, but it should be much more continuous. 

To partly account for that we have added, in Fig. |4] and 
in Fig. a constant luminosity equal to 0.1 per cent of the 
Eddington value. The thick solid lines are the sum of this con- 
tribution and the contribution produced by the shell-shell col- 
lisions. 

The three different curves in Fig. 0] show the effect to 
change the time interval between the launching of blobs. For 
increasing time intervals, less spikes per unit time are produced 
and the overall process becomes less efficient, since there are 
less collisions. 

Note that, in agreement with our simple estimates, the av- 
erage emitted luminosity (at least for short At) is around a few 
per cent of the jet initial kinetic power. 

Fig-Elshows the optical depth tt, the electron temperature 
kT and the Comptonization y parameter as a function of the lu- 
minosity produced in each collision (i.e. this is not the sum over 
the "on" shells used in Fig. |4j. For reference, we have drawn 



a vertical line corresponding to 0. 1 per cent of the Eddington 
luminosity and corresponding to the stationary component as- 
sumed above. For this case we have assumed At = 4i? g /c as 
in the top panel of Fig.|4] As can be seen, when the luminosity 
is relevant (i.e. L° bs ~ 10 _1 Lj = 10 -2 LEddX we nave optical 
depths between 0.5 and 2, temperatures around 100 keV and 
y ~ 1 . For this case we obtain values of the optical depth which 
are a few times larger than the initial ones, due to the previous 
collision done by the shell, as shown in Fig.|6] 

5.1. Observed spectra 

Since we can calculate the optical depth and temperature for 
each shell-shell collision, we can calculate the emitted spec- 
trum assuming thermal Comptonization as the radiative pro- 
cess. To this end we use the analytical formulae of Titarchuk & 
Mastichiadis (1994), and assume that the soft radiation field is 
a blackbody peaking at some energy eo- F° r all cases discussed 
in the following, we have assumed eo = 5 eV. For each observ- 
ing time, we sum up the spectra of the "on" shells, including 
the contribution of the steady component, whose bolometric 
luminosity amounts to the 0. 1 per cent of the Eddington lumi- 
nosity. We have assumed that the spectral shape of this steady 
component is F(v) oc v~ 1 e~ v ^ c , with hv c = 150 keV. 

This specific chosen spectral index may account for the 
contribution of the hot corona, above the accretion disk, while 
it is not clear if it can be directly associated to the jet genera- 
tion process (e.g., to a particle injection in the jet smoother than 
the assumed on-off mechanism). In this latter case the usual 
feedback process operating in the hot corona model, fixing the 
spectral index close to unity, does not work. In fact, when the 
power of the steady component is much smaller than the lumi- 
nosity of the disk, there is no feedback, and the Comptonization 
spectrum should be much steeper than unity. 

Fig- |U shows some example of spectra born out from our 
simulations. They correspond to the case shown in the top 
panel of Fig. |4] assuming a black hole of 10 s solar masses. 
Spectra are calculated every 40 seconds, but are shown ev- 
ery 800 seconds (corresponding to 0.53i? g /c), for clarity, from 
t/(R g /c) = 240 to t/(R g /c) = 300. Fig.|5]shows again the 
light curve corresponding to the top panel of Fig. @] together 
with the associated 1-10 keV spectral index, d[ 1 _ 10 ] and the 
light curve of the 1-10 keV flux. Fig.llOlshows the histograms 
of a[i_io] and the softer spectral index a[o.i-i] (between 0.1 
and 1 keV). 

By construction, a spectral index different from unity cor- 
responds to the jet emission dominating over the steady com- 
ponent. Most of the time, the observed slopes corresponding to 
the jet emission are steeper than unity (as shown in Fig. 1 1 Oi . 
and even more so at lower frequencies. 

5.2. Trends 

In this section we summarize the results of the survey of the 
parameter space aimed to test the sensitivity of our results to 
our assumptions. We indicate the values of the parameters we 
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log E [keV] 

Fig. 8. Example of time dependent spectra for the case illus- 
trated in the top panel of Fig.|4] assuming a black hole of 10 8 
solar masses. Spectra are calculated every 2.7 x 10~ 2 i? g /c, but 
for clarity the figure shows spectra every 0.53i? g /c (i.e. ev- 
ery 800 seconds for a black hole of 10 8 solar masses), from 
t/(R g /c) = 240 and 300. Vertical line mark the 1-10 keV 
band. 



have changed, all other parameters have values equal to the 
case shown in the top panel of Fig.|4] 

- Initial velocity — Fig. shows the effect of changing the 
average initial velocity of the ejected shell ((/3o) =0.3, 0.5 
and 0.7 from top to bottom). It can be seen that an increase 
of the average initial velocity has the effect to increase the 
density of spikes and their average luminosity. Increasing 
the average velocity, in fact, increase the average number of 
collisions per unit time and the average efficiency of each 
collision. The overall process becomes more efficient. 

- Launching site — Increasing the launching site from the 
black hole (we tried Rq = 3, 5 and 7 i? g ) the bolometric 
and the 1-10 keV light curves becomes more spiky, with 
the jet dominant (in the 1-10 keV band) an increasing num- 
ber of times over the steady component. The net effect is, 
basically, the same as increasing the initial velocity. 

- Shell size — Increasing Arj (we tried Arj = 1,2 and 4 i? g , 
always with rj — 4i? g ), the light curve becomes slightly 
less spiky, as the duration of the collisions is proportional 
to Arj. 

- Spin of the black hole — We have found negligible 
differences between a maximally rotating Kerr and a 
Schwarzchild black hole for the launching site Rq = 5i? g . 
There is still a negligible difference for the ca.se Rq — 3_/?g 
and/? = 0.7. 

- Ratio of jet to disk power — We consider this as the 
most important parameter. As the jet power increases, the 
jet emission becomes (obviously) more and more dominant 
over the steady emission, and the average spectral index 




100 200 300 400 

Time [R G /c] 



500 



Fig. 9. Light curves of the integrated X-ray luminosity (in 
units of the Eddington one) corresponding the top panel of Fig. 
[4] the light curve of the spectral index and flux between 1 and 
10 keV. We have assumed that the spectral shape of the steady 
component is F(y) oc v~ x e~ v ' Vc . 




Fig. 10. Histograms of the soft (0.1-1 keV) and hard (1-10 
keV) X-ray spectral indices for the case shown in Fig. [9] 



flattens, to become flatter than unity for Lj ~ 20La. For 
larger values of the jet power, the radiation produced by the 
jet is larger that the luminosity radiated by the disk. The ra- 
diation coming from the jet and reprocessed by the disk can 
no longer be neglected, and our treatment becomes invalid. 
However, in this case, the power in the optical-UV and in 
the X-ray bands are comparable, in contrast with the ma- 
jority (albeit not all) of Seyfert galaxies. 
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6. Summary and discussion 

We have explored the possibility that all Active Galactic Nuclei 
form jets or outflows, with a range of velocities, but with a 
power which is comparable to the power extracted by the ac- 
cretion process. We have then assumed that in most cases (cor- 
responding to radio-quiet sources) the jets are launched with 
velocities smaller than the escape velocity. Mainly for simplic- 
ity, but also in analogy with the "internal shock" scenario pro- 
posed to work in radio-loud sources and in gamma-ray bursts, 
we have further assumed that the central engine works inter- 
mittently, producing shells or blobs. A shell with (3 < f3 csc 
will reach a maximum distance from the hole, then stop and 
invert its motion, and may eventually collide with the succes- 
sive shell. In this case the bulk kinetic energy of the two shells 
is dissipated, and the fraction of it which is given to electrons 
can be transformed into radiation. The accelerated electrons are 
embedded in the dense radiation field produced by the accre- 
tion disk: they cool rapidly by the inverse Compton process, 
producing the X-ray continuum. Simple energy balance is suf- 
ficient to estimate the Comptonization parameter as a function 
of the power dissipated by the colliding blobs and the disk lu- 
minosity. By construction, the shells are moving with velocities 
smaller than the escape speed, yet they carry a power which 
is comparable to that extracted by accretion. This implies that 
the shells are "heavy". They cannot be formed by electron- 
positron pairs only, because in this case the corresponding op- 
tical depth is so large that most of them annihilate in less than 
a dynamical time. The required density in protons and the ac- 
companying electrons is large enough to limit the importance 
of pairs not only as energy carriers, but also as scatterers. 

The main aim of this paper is to investigate the general 
properties of the proposed idea, to check if it can work at least 
at the first order of approximation. 

Our "aborted jet" scenario is not necessarily alternative to 
the popular "disk-corona" model. Both processes could be ac- 
tive and contribute to the formation of the high energy contin- 
uum in the same source. On the other hand we would like to 
stress that in our proposed scenario the source of energy could 
be the spin of the hole, besides accretion. Pushing this possi- 
bility to the limit (i.e. all the high energy emission produced by 
AGNs comes from the rotational energy of their black hole), 
would result in the remarkable fact that it is the black hole 
spin, rather than accretion, which produces the bulk of the X- 
ray background. It is then instructive to isolate the "aborted jet" 
process in order to find ways to confirm or falsify this scenario. 

One of the clearest difference with the disk-corona model 
is that the dissipation of energy should occur along the axis 
of rotation of the black hole. This implies that the X-ray flux 
coming from the colliding shells will illuminate preferentially 
the inner part of the disk, especially when they collide close 
to the hole. This may solve the problem of the formation of 
the strong red wings of the relativistic iron line observed in 
MCG-6-30-15 (Wilms et al. 2001), which requires an "illumi- 
nator" emissivity strongly increasing towards the black hole. It 
should be noted that we have neglected, for simplicity, the light 
bending due to the strong gravity (see e.g. Martocchia et al. 
2002), which results in an enhanced illumination of the inner- 



most disk regions. The illumination could be further enhanced 
by anisotropic Compton scattering (since the seed photons are 
coming from the disk, more inverse Compton radiation is chan- 
neled back towards the disk than along the viewing angle, see 
e.g. Ghisellini et al. 1991; Malzac et al. 1998). Another cause 
of anisotropy is beaming of the X-ray radiation, which is pref- 
erentially emitted towards the accretion disk in efficient shell- 
shell collisions. According to our simulations, in fact, the most 
efficient collisions are between massive blobs coming back to 
the disk and having already experienced some collisions, and 
newly generated blobs moving in the opposite direction. While 
this may help explaining the large equivalent widths of Fe lines 
observed in a few cases (notably MCG-6-30-15), it is appar- 
ently at odds with the relative paucity of relativistic iron lines 
observed by XMM-Newton (e.g. Reeves et al. 2003). However, 
it should be noted that the increase of illuminating X-ray pho- 
tons may result in a significant ionization of the innermost re- 
gions of the accretion disc, making predictions on the iron line 
intensity less straightforward (e.g. Nayakshin & Kazanas 2002, 
and references therein). Detailed calculations of the iron line 
properties are beyond the scope of this paper, and are deferred 
to a future work. 

We note that an important piece of information may come 
from the observations of the Compton reflection continuum and 
iron line in radio-loud sources. If the jet in these sources is 
successful, in fact, it should not illuminate much of the ac- 
cretion disk, and therefore these objects should have weaker 
reflection features produced by the corona only. This seems 
indeed to be the case (e.g. Grandi et al. 2002, and references 
therein). Then the equivalent width of the fluorescent iron lines 
in radio-galaxies may measure the importance of the corona 
with respect to the jet in producing the thermal X-ray contin- 
uum, once the data are purified from all other additional con- 
tributions (e.g. the non thermal radiation from the jet). 

The spectral index of the jet emission, calculated in our 
simulations, is generally steeper than the average spectral in- 
dex observed in Seyfert galaxies (i.e. a x ~ 1). As it is, our 
model requires therefore the presence of a steadier component, 
with the "right" spectral index, contributing to X-ray band. 
This steady component should have a bolometric luminosity 
which is, on average, smaller than the average power of the 
jet, even if its relative contribution in the 1-10 keV band is 
more important. The steep jet emission, when contributing no- 
tably to the 1-10 keV band, would steepen the overall spectral 
index and increase the flux. It would then produce a "steeper 
when brighter" behavior as observed in Seyfert galaxies (e.g. 
Zdziarski et al. 2003). Occasionaly, instead, the jet emission 
is both dominant and characterized by a fiat spectrum, and we 
have then the opposite behavior, i.e. "harder when brighter", 
but this occurs more rarely. 

As it is, our model explains the X-ray properties of Narrow 
Line Seyfert 1 galaxies (Boiler et al. 1996; Brandt et al. 1997; 
Cancelliere & Comastri 2002). These sources are in fact char- 
acterized, on one hand, by a 2-10 keV spectral index between 1 
and 1 .5 (and an even steeper spectrum in the softer band), and, 
on the other hand, by a short term, large amplitude variability. It 
is then possible than the main difference between Narrow Line 
Seyl (including in this class also sources like MCG-6-30-15 
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which have broad lines but in X-rays behaves like NLSyls) 
and classical Seyferts is the ratio between jet and disc/corona 
emission. It is worth noting that NLSyl are widely believed 
to have a larger L/L^dd ratio than classical Seyferts, which 
again can be explained by an enhanced jet emission. If this is 
true, one could speculate that the physical parameter behind 
the NLSyl X-ray behaviour is not (or at least not only) the ac- 
cretion rate, as usually supposed, but the presence of a more 
powerful aborted jet. 

Regarding broad line, classical Seyfert 1 galaxies, we 
should however consider that our model neglects, for simplic- 
ity and ease of calculation, a few important physical effects. 
One of these concerns light bending, important when the emit- 
ting spot is very close to the black hole. This effect is expected 
to change the observed X-ray luminosities only by a factor 
of a few (Martocchia et al. 2002), but that can nevertheless 
be very important for a detailed study of spectral evolution. 
Indeed, a different degree of light bending and gravitational 
redshift corresponding to different heights of the illuminator 
above the black hole can explain the puzzling temporal behav- 
ior of MCG-6-30-15, where the continuum and iron line vari- 
abilities are decoupled (Miniutti et al. 2003). 

Another important effect, neglected here, is the feedback 
between the luminosity produced by the jet and the disk emis- 
sion, important for large ratios between the jet and the disk 
powers. In these cases the radiation reprocessed by the disk 
can become important and introduce the same kind of feedback 
which makes the hot corona model to work, producing spectral 
indices close to unity in the X-ray band. 

Finally, we would like to comment about the difference be- 
tween radio-loud and radio-quiet sources. In our scenario, this 
is mainly a difference in mass loading, coupled with a possi- 
ble difference in jet power. The central engine in radio-loud 
sources succeeds in accelerating jets at speeds larger than the 
escape velocity: in these sources the jet power can dominate the 
total energetics (as in BL Lac objects), and the outflow mass 
rate is of the order of a per cent of the accretion rate. The "jet" 
of radio-quiet sources may not be much less powerful than in 
radio-loud objects, if it contributes significantly to the forma- 
tion of the X-ray flux. What should be different is the outflow- 
ing mass rate, which must be greater in radio-quiet objects, 
making their "jets" move slower. If all jets are powered by the 
extraction of rotational energy from a spinning black hole, it is 
then possible that it is this mechanism, and not accretion, to be 
responsible for all the high energy radiation produced by AGNs 
(i.e. all the X-ray and the 7-ray flux). 

It is also possible that a specific source, usually radio- 
quiet, occasionally may launch "successful" shells, with rel- 
ativistic speeds. However, these "successful jet episodes" must 
be rare in AGN, since we rarely see "fossil" long lived weak 
radio lobes in not jetted sources. This may occur more of- 
ten in galactic micro-quasars, and be associated with the ma- 
jor radio-flares. The bulk Lorentz factor associated with ma- 
jor radio events in GRS 1915+105 is relatively small, perhaps 
suggesting that, when radio-weak, the jet is not successfully 
launched because it does not attain bulk speeds larger than the 
escape velocity. These sources may therefore be the "missing 



link" between radio-loud and radio-weak sources, changing 
from time to time from one class to the other. 
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